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Abstract 

Minimum-time accelerations of the F100 turbofan 
engine are presented. A piecewise-linear engine 
model, having three state variables and four con- 
trol variables, is used to obtain the minimum-time 
solutions. The linear model which applies at a 
given time in the trajectory is determined by cal- 
culating a normalized "distance" from the current 
state to the equilibrium state associated with each 
linear model. The linear model associated with the 
closest equilibrium point is then used. The con- 
trol histories for the minimum-time solutions are 
used as input to a nonlinear simulation of the F1Q0 
engine to verify the accuracy of the piecewise- 
linear solutions. 


INTRODUCTION 

Modern, high-performance turbojet and turbofan 
engines are generally equipped with one or more 
variable geometry features in order to provide max- 
imum propulsive efficiency over a range of engine 
power settings and flight conditions. For example, 
the J-85 engine (a one-spool turbojet used in the 
F5 aircraft) has variable inlet guide vanes, and 
variable bleeds in three stages of its eight-stage 
compressor. The TF30 engine (a two-spool turbofan, 
used in the Fill and F14 aircraft) has variable 
bleeds in the low and high compressors. The F100 
engine (a two-spool turbofan, used in the F15 and 
F16 aircraft) has variable fan inlet guide vanes 
and variable compressor stator vanes. Each of 
these engines also has a variable-area exhaust noz- 
zle and an afterburner. Variable area turbines, 
although not yet in operational use, have been 
tested on technology demonstrator engines. 

Propulsive efficiency is probaoly the most im- 
portant measure of an aircraft engine's performance. 
However, another important measure is the time re- 
quired to accelerate from one thrust level to 
another. Engine acceleration is one of the func- 
tions of the engine control system, and may be ac- 
complished via open-loop scheduling or closed-loop 
control. For each of the three engines referred to 
above, engine acceleration is accomplished by con- 
trolling fuel flow. The variable geometry features 
are not utilized in this function. Instead, these 
variables are scheduled in a steady-state manner to 
maximize propulsive efficiency. 

In recent years, linear-quadratic regulator 


theory has been developed for the design of multi- 
input multi-output control systems. An account of 
the theory and application is given, for example, 
in reference 1. Use of the theory has been facili- 
tated by computer programs such as those de- 
scribed in references 2 and 3, which rapidly and 
efficiently calculate the optimal feedback control 
gains, given the system description and perform- 
ance index. This theory has been applied recently 
to the design of control systems for aircraft gas 
turbine engines. In addition to the design of 
regulators, the problem of minimizing acceleration 
time has also been considered. 

Michael and Farrar (refs. 4 and 5) applied 
linear quadratic regulator theory to the design of 
controls for the F401 turbofan engine. The non- 
linear system equations were linearized about five 
different equilibrium points, a quadratic perform- 
ance index intended to minimize acceleration time 
was formulated, and feedback control gains were 
determined. A nonlinear feedback control law was 
developed by curvefitting the resulting control 
gains as a function of compressor speed. 

Weinberg (ref. 6) applied linear-quadratic 
regulator theory to the design of controls for the 
F100 engine. He showed that this engine could be 
adequately represented by three state variables - 
fan speed, compressor speed, and augmentor pres- 
sure, Four control variables were utilized, and 
linearized engine models were obtained at two 
equilibrium points. The problem of minimizing ac- 
celeration time was considered, and control system 
gains were derived by conducting small perturba- 
tion optimizations at each of two equilibrium 
points, using a quadratic performance index. The 
control gains were switched at a fixed value of 
fan speed, rather than varied in a continuous man- 
ner as in references 4 and 5. 

In reference 7, Sevich and Beattie considered 
the minimization of acceleration time for a turbo- 
jet engine, using fuel flow and exhaust nozzle 
area as control variables. They used a quadratic 
performance index to approximate minimum-time solu- 
tions, as in references 4 to 6. However, they used 
a nonlinear engine model, rather than a series of 
linear models. The result was an open-loop, opti- 
mal trajectory. The controls were assumed to be 
piecewise constant, and the performance index was 
minimized by using a conjugate-gradient search 
technique. 

Dehoff et al. (ref. 8) used linear-quadratic 
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regulator theory to design controls for the F100 
engine. The control gains were generated using 
linear models with five state variables and four 
control variables at several equilibrium points. 
Principal emphasis was on the regulator design. 
Although acceleration control was considered, there 
was no specific attempt at minimizing the accelera- 
tion time. 

References A through 8 all made use of inte- 
gral, quadratic performance indices, in which both 
state and control deviations from some desired tra- 
jectory were penalized. The coefficients of the 
penalty terms were adjusted in an attempt to mini- 
mize the acceleration time. However, none of these 
reports can claim their final histories produce 
truly minimum-time accelerations. 

In this paper, minimum-time acceleration his- 
tories are presented for the F100 turbofan engine. 
Four control variables, i.e., fuel flow, exhaust 
nozzle area, inlet guide vane position, ar.d com- 
pressor stator vane position, are utilized. A 
piecewise-linear engine model having three state 
variables and four control variables was used to 
obtain the minimum-time solutions (ref. 9). The 
linear models were obtained by linearization of a 
nonlinear model at four sea-level, static equilib- 
rium points, and were taken from reference 10. The 
linear model which applied at a given time in the 
tr« jectory was determined by calculating a normal- 
ized "distance” from the current state to the state 
at each of the equilibrium points; the linear model 
associated with the closest equilibrium point was 
then used. Linear state/control constraints which 
correspond to speed, temperature, pressure, and 
control limits were considered. The minimum-time 
solutions are used as inputs to a nonlinear com- 
puter simulation of the F100 engine (ref. 11) to 
verify the accuracy of the piecewise-linear solu- 
tions. 

The minimum-time solutions were obtained by 
using a new algorithm, which is described in refer- 
ences 9 and 12. A more complete discussion and 
presentation of the results, as well as the rele- 
vant mathematical formulation, may be found in ref- 
erence 9. 

ENGINE MODEL 

In a recent contractual effort under joint 
Air Force/NASA sponsorship (ref. 8), Systems Con- 
trol Inc. (SCI) used linear quadratic regulator 
theory to design controls for the F100 engine. 
Linear models having 16 states and reduced-order 
models having five states were provided to SCI by 
Pratt & Whitney Aircraft (P&WA) for a number of 
equilibrium points at different flight conditions 
and power settings (PLA) . Some of these models are 
given in reference 10. Four reduced-order, sea- 
level static models from the P&WA/SCI study, 

(PLA - 36, 52, 67, and 83 degrees) were used to 
generate the results presented herein. 

In order to determine if the fifth-order 
models could be further reduced, the eigenvalues 
and eigenvectors of the fifth-order models were 
calculated. It was found that there were three 


dominant eigenvalues , one real and one complex 
pair. The other two eigenvalues had much larger 
values (higher frequency). Therefore, the number 
of states was reduced from five to three, as de- 
scribed in reference 9. 

TRANSIENT PERFORMANCE 

The problem of minimizing the time required to 
accelerate the F100 engine from one equilibrium 
thrust level to another is considered. In solving 
this problem, a piecewise-linear model of the F100 
engine, operating at the sea-level, static condi- 
tion, was constructed from the four equilibrium 
linear models at PLA - 36, 52, 67, and 83 degrees. 
At any given time, the model whose equilibrium 
state is "closest" to the actual state at that 
time is used to represent the engine. 

During an acceleration from near-idle to 
intermediate thrust, the engine is first repre- 
sented by the PLA * 36° model. As the engine 
accelerates, the model switches successively to 
the 52, 67, and 83 degree models. 

There are path inequality constraints which 
must be satisfied along a trajectory. Some of 
these constraints correspond to engine physical 
limits, others to control mechanical limits. The 
following constraints are assumed for this paper. 

(1) Turbine inlet temperature cannot exceed 
the equilibrium value at intermediate thrust by 
more than 50 degrees R. 

(2) Fan and compressor speeds cannot exceed 
the equilibrium values at intermediate thrust by 
more than 50 rpm. 

(3) Fan and compressor surge margins must not 
be less than 5 percent. 

(A) Inlet guide vanes, compressor vanes, ex- 
haust nozzle area and fuel flow rate must not ex- 
ceed their limits. 

STATEMENT OF THE PROBLEM 

The problem to be solved is to find the values 
of the four control variables, as functions of 
timn , which minimize the "terminal error" (dis- 
tance from the final state to the desired final 
state) for a specified acceleration time. Tills 
must be accomplished while satisfying the system 
equations and path constraints. A sequence of 
solutions to such problems may be used to find the 
minimum-time solution for a given value of term- 
inal error. The mathematical problem formulation, 
necessary conditions for optimality and solution 
procedure may be found in reference 9. 

RESULTS 

The problem of minimizing the terminal error 
for an acceleration from near-idle (PLA - 2A°) to 
intermediate thrust in 0.75 second was considered 
in reference 9. The responses of the problem var- 
iables (states, outputs, and controls) for the 
optimal solution are shown in figure 1. The state 
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variables (fan speed, compressor speed, and aug- 
mentor pressure) are shown in figures 1(a) to (c) , 
respectively. It can be seen that the states 
approach the desired final values smoothly and 
with no overshoot. 


COMPARISON OF NONLINEAR AND 
PIECEWISE- LINEAR RESULTS 


The output variables are shown in figures 1(d) 
to (J). Because a piecewise-linear model has been 
utilized, the outputs are in general discontinuous 
at model switching points as well as at points of 
discontinuous control. If the piecewise linear 
model is a good representation of the engine, the 
discontinuities in the outputs at model switching 
points should be small. In figure 1, the disconti- 
nuities in thrust, low-pressure turbine inlet tem- 
perature, and airflow are so small that they are 
not visible. However, the discontinuities are sub- 
stantial for fan and compressor surge margins and 
for combustor pressure. 

The optimal control strategy results in the 
high-pressure turbine inlet temperature having its 
maximum value for the entire trajectory; this is 
shown in figure 1(d). Fan and compressor surge 
margins (figs. 1(e) and (f)) remain well above ac- 
ceptable mlnimums. Low-pressure turbine inlet tem- 
perature (fig. 1(g)) is very nearly constant. 

Thrust (fig. 1(h)) increases smoothly and mono- 
tonlcally. 

The optimal control histories are shown in 
figures l(k) to (n) . Fuel flow jumps initially 
from its idle value, then increases slowly. The 
optimal value of nozzle area is constant, and equal 
to its lower mechanical limit. For the accelera- 
tion shown in figure 1, the values of inlet guide 
vane position (IGV) and compressor variable vane 
position (HVS) were constrained to be within ±7 de- 
grees of the bill-of -material (BOM) control sched- 
uled values. This was necessary because of limited 
off-schedule test data and to avoid violating 
flutter boundaries. For this case, IGV switched 
from 7 degrees below the scheduled value to 7 de- 
grees above the scheduled value at 0.56 second; the 
value of HVS was 7 degrees below the schedule for 
the entire trajectory. 

Figure 2 shows the terminal error as a func- 
tion of acceleration time lor accelerations from 
near-idle (PLA - 24°) to intermediate thrust 
(PLA ■ 83°). A terminal error of 0.10, for exam- 
ple, indicates that all three states are within 
10 percent of their desired final values at the 
completion of the trajectory. Results are pre- 
sented for values of IGV and HVS which are fully 
optimized, scheduled, and scheduled ±7°. Figure 2 
shows, for example, that an acceleration time of 
0.80 second is required to reduce the terminal 
error to 0.05 when scheduled IGV and HVS are used. 
However, when the IGV and HVS are controlled opti- 
mally within ±7 degrees of the scheduled values, 
the Acceleration time is reduced to 0.74 second. 
Fully optimized IGV and HVS resulted in an acceler- 
ation time of 0.65 second. The corresponding tra- 
jectories have the same characteristics as shown in 
figure 1. In particular, the high-pressure turbine 
inlet temperature has its maximum value throughout 
each of the trajectories. 


In order to determine the accuracy of the 
piecewise-linear solutions, the control histories 
for the time optimal solutions were used as input 
to a nonlinear F100 engine simulation (ref. 11). 
Figure 3 presents a comparison of nonlinear and 
piecewise-linear transient responses. The control 
variable histories are based on an optimal accel- 
eration from PLA ■ 24 degrees to intermediate 
thrust, with a specified acceleration time of 
0.8 second. For this case, the values of IGV and 
HVS were constrained to be within ±7 degrees of 
the scheduled values. It can be seen that the 
nonlinear and piecewise-linear responses of com- 
pressor speed and augmentor pressure are in good 
agreement. However, differences are observed in 
the fan speed responses. Also, there are substan- 
tial and Important differences in the high- 
pressure turbine inlet temperature and fan and 
compressor surge margin responses. The nonlinear 
results show that the maximum value of the high- 
pressure turbine inlet temperature is violated by 
a large amount , and that the compressor surges at 
about 0.06 second. The fan does not surge, but 
the fan surge margin does fall below the minimum 
value of 5 percent early in the trajectory. 

There are several possible explanations for 
the differences between nonlinear and piecewise- 
linear results observed in figure 3. They are as 
follows : 

(1) The individual linear models may not be 
precisely correct even for small perturbation in- 
puts. 

(2) The number of equilibrium linear models 
may be inadequate to accurately describe the sys- 
tem nonllnearltles with respect to the state var- 
iables. 

(3) There may be substantial nonlinearities 
with respect to the control variables, which are 
not Included in the piecewise-linear model. 

(4) Model reduction from sixteenth order to 
third order may have resulted in modeling inaccu- 
racies. 

An investigation is currently underway to deter- 
mine and eliminate the cause of the piecewise- 
linear model inaccuracies. 

COMPARISON OF MINIMUM-TIME AND 
BOM CONTROL RESPONSES 

It is of interest to compare the transient re- 
sponses obtained by using optimal, minimum-time 
strategy with those obtained by using the BOM con- 
trol. Such a comparison is made in figure 4, for 
an acceleration from PLA ■ 24 degrees to inter- 
mediate thrust. The minimum time controls were 
obtained by using the piecewise-linear model with 
the values of IGV and HVS constrained to be within 
±7 degrees of the BOM scheduled values. However, 
the results presented in figure 4 were obtained by 
using the minimum-time and BOM controls in con- 
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Junction with the nonlinear F1Q0 engine simulation. 
The results show that the minimum-time control 
strategy produced a more rapid acceleration to 
Intermediate thrust. It appears that the principal 
reason for the improved acceleration is the much 
more rapid Increase in fuel flow, which results in 
a rapid increase in high-pressure turbine lnlit 
temperature. Naturally, the comparison of perform- 
ance is invalidated because of the violation of 
constraints which occurs. Nevertheless, it seems 
highly probable that substantial improvement in ac- 
celeration time can be made without violation of 
engine constraints since the high-pressure turbine 
inlet temperature (fig. 4(d)) Increases very slowly 
when the BOM control is used. 

CONCLUDING REMARKS 

In this paper, minimum-time acceleration his- 
tories are presented for the F100 engine. A 
piecewise-linear engine model, having three state 
variables and four control variables, was used to 
obtain the minimum time solutions. The resulting 
control histories are used as Inputs to a nonlinear 
simulation of the F100 engine to verify the accu- 
racy of the piecewise-linear solutions. 

A comparison of the nonlinear and piecewise- 
linear solutions revealed significant differences 
in some of the engine responses. Several possible 
explanations for these differences were noted; this 
problem is currently under investigation. 

In this pLper, it was assumed that the control 
variables could jump dlscontinuously from one value 
to another. Studies are currently underway to find 
minimum-time acceleration histories which satisfy 
the control variable rate limits in addition to the 
other constraints. 

Results presented herein indicate that im- 
proved steady-state and transient performance may 
be obtained by using optimal control strategy. 

Such strategy sometimes calls for operation of the 
controls in a manner which has not been previously 
tested experimentally, even though the nonlinear 
simulation contains equations for predicting the 
effect of such control action. Further experi- 
mental testing is indicated in order to systemati- 
cally explore the steady-state and transient ef- 
fects of off-scheduled control action, and to de- 
termine if the predicted performance gains can be 
achieved. 
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